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ABSTRACT 



The search for gamma-ray line emission from the ra- 
dioactive decay of 60 Fe figures among the prime sci- 
entific objectives of the INTEGRAL mission. 60 Fe is 
believed to be primarily produced in core-collapse su- 
pernovae, although other sites, such as carbon defla- 
gration supernovae or intermediate mass AGB stars 
have also been suggested. We present first results 
of our search for the 1173 and 1332 keV gamma-ray 
lines of 60 Fe decay in the SPI data of the first mis- 
sion year. So far we can only report upper flux limits 
whose levels are determined by systematic uncertain- 
ties in the treatment of the instrumental background. 
Once we understand the background uncertainties 
better, the SPI data should provide stringent con- 
straints on 60 Fe nucleosynthesis in our Galaxy. 



Key words: gamma rays: observations, lines; super- 
novae; nucleosynthesis. 



1. INTRODUCTION 



Searching for the gamma-ray lines at 1173.23 and 
1332.49 keV that arise from the radioactive de- 
cay chain 60 Fe — * 60 Co — * 60 Ni figures among 
the prime objectives of the INTEGRAL mission. 
60 Fe is believed to be primarily produced in core- 
collapse supernovae, although other sites, such as 
carbon deflagration supernovae or intermediate mass 
AGB stars have also be en suggested (see review in 
iKnoHlseder fe Vedrermel. 2001). 

Theoretical estimates of 60 Fe yields from core- 
collapse supernovae suffer from considerable uncer- 
tainties in nuclear reactions and in sup ernova ex- 
plosion model details which affect these (|Prantzost 



l2004j) . hence observing the 60 Fe gamma-ray lines will 
provide valuable constraints on supernova nucleosyn- 
thesis. From comparisons of isotope abundance mea- 
surements whic h are co-produced with 60 Fe, it has 
been estimated iLeising fe Share) . I1994D that a total 
amount of 0.9Mq should exist in the present Galaxy, 
leading to a flux of ~ 3 x 10~ 5 ph cm _2 s _1 in each 
of the lines (corresponding to about 8% of the flux 
in the 1809 keV line that arises from the radioactive 
decay chain 26 A1 -> 26 Mg). Recently. ISmithl f2004b) 
has reported a first detection of the 60 Fe lines by 
the solar observatory RHESSI, yet during this con- 
ference he revised the si gnificance of th e result to 
a marginal level of 2.6a llSmithl . l2004a|) . By com- 
bining both lines, Smith derives a flux in each line of 
(5.2±2.0) x 10 -5 ph cm~ 2 s~ 1 from the inner Galaxy, 
corresponding to about 10% of the flux that RHESSI 
observes in the 1809 keV line. 

In this paper we present our first results from the 
search for 60 Fe decay emission in the data of the spec- 
trometer SPI aboard INTEGRAL. For the moment 
we have no indications of significant 60 Fe line emis- 
sion, yet we are still facing considerable problems in 
the treatment of the instrumental background that 
have to be resolved before SPI can reach its actual 
sensitivity limit. 



2. ANALYSED DATA 



The data that we analysed in this work consist 
of the INTEGRAL first year's spring and autumn 
Galactic Centre Deep Exposures (GCDE1 and 2), 
with 1.8 and 2.0 Ms of exposure, respectively, cov- 
ering galactic longitudes I = ±30° and latitudes 
b = ±20°. GCDE1 includes INTEGRAL orbits 47- 
66, GCDE2 includes orbits 97-123. Only the analysis 
of single-detector event data is presented in this pa- 
per. The analysis of multiple-detector event data has 
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Figure 1. SPI instrumental background spectrum (crosses) in the area around the 60 Fe line energies together 
with the background model (histogram). The summed and degradation corrected data for GCDE1 and 2 are 
shown. The grey areas (red in the colour version of the paper) indicate the energy bands that have been used 
for modelling the continuum emission underlying the line (indicated by the light line). Instrumental background 
lines have been labelled by their source isotopes, where identified. Note that we only attempt to model the 
60 Co background lines, while other background lines which obey different time variabilities are not described 
satisfactorily. 



been started, yet due to issues concerning the back- 
ground modelling we have so far focused our atten- 
tion on the single-detector event data. Energy cali- 
bration has been performed orbit-by-orbit, resulting 
in a relative (orbit-to-orbit) calibration precision of 
~ 0.01 keV and an absolute accuracy of ~ 0.05 keV 
<)Loniou et alll2004|) . 



3. BACKGROUND MODELLING 



Background modelling is performed in our usual ap- 
proach, where a two component model aims to de- 
scribe both the instrumental line and continuum 
components. At the rest energies of the astrophysical 
60 Fe lines (1173.23 and 1332.49 keV) the SPI back- 
ground (cf. Fig. ^) shows two instrumental line fea- 
tures that have been identified as activation lines of 



the radio-isotope 60 Co (i.e. t he same isotope that we 
aim also to observe in space; IWeidenspointner et all 
2003). 

For the line component, the spectral shape and the 
relative counting rates in each detector have been 
estimated using empty field observations, assuming 
that both do not vary in time (in fact, detector degra- 
dation leads to a modification of the instrumental 
line shape s that is corrected for after background 
modelling; iKnodlseder et al.1 2004). A total empty 
field exposure of 4.1 Ms has been collected by gath- 
ering all data for which SPI pointed at galactic lati- 
tudes \b\ > 20° (covering INTEGRAL orbits 14-106). 
To extract the line component, a linear continuum 
background has been subtracted from the spectra, 
estimated from adjacent energy bands (1163-1169 & 
1177-1183 keV for the continuum under the 1173 keV 
line and 1318-1329 & 1351-1360 keV for the contin- 
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uum under the 1332 keV line). 

The time variability of the line component has been 
modelled by what we call a background variation 
template. This template is a combination of back- 
ground activity tracers that have been adjusted to 
explain the line background variation in the empty 
field observations. The idea behind this approach 
is that an instrumental background line has in gen- 
eral N physical contributions which may be prompt 
(instantaneous) or may eventually come from decay 
chains with various decay constants r. Thus, in 
general, the counting rate of the instrumental back- 
ground line component can be written as 

where the index i runs over the N contributions to 
the background line, ri(t) is the counting rate of the 
activity tracer used to model the contribution i, tj is 
the decay constant of contribution i, and a, is the rel- 
ative amplitude of contribution i. The time integral 
is taken from the start of the INTEGRAL mission 
up to the date of interest t. 

For the 60 Fe analysis, we tried different combinations 
of activity tracers such as the germanium saturated 
event rate as a tracer of the high-energy particle irra- 
diation of the SPI telescope or the plastic scintillator 
event rate as a tracer of the low-energy particle irra- 
diation. As example, we show in Fig. a template 
made of two contributions, which consist of a prompt 
background component (i.e. r = 0), modelled using 
the germanium saturated event rate, and a 60 Co ac- 
tivation component, modelled using the germanium 
saturated event rate convolved with an exponential 
law with a time constant of t = 7.6 yr. The rel- 
ative amplitudes at of each contribution have been 
adjusted using the empty field data. Replacing the 
germanium saturated rate by the plastic scintillator 
event rate or adding further components does not 
alter our background model significantly. Therefore 
we only show results in this paper that have been 
obtained by using the template of Fig. [3 
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Figure 2. Background variation template of the in- 
strumental background line at 1173.23 keV. 



The continuum background has been estimated by 
adjusting linear functions to the data in energy bands 
adjacent to the line (the same energy bands that 
have been used for the line component extraction in 
the empty field data). In principle, this adjustment 
should be done for each detector and pointing, yet 
the small number of counts in the dataspace requires 
an adaptive rebinning (in the time axis), request- 
ing a minimum number of 1000 events per adjacent 
energy bin. This rebinning has been found to be cru- 
cial in reducing systematic uncertainties in the line 
shape analysis. As long as the number of events used 
to determine the continuum level is sufficiently large 
(> 1000), the precise smoothing parameters have lit- 
tle impact on the result. 

Figure shows the background subtracted spectrum 
for the summed GCDE1 and 2 data. As we at- 
tempted only to model the instrumental 60 Co lines, 
strong residuals remain at the energies of instrumen- 
tal background lines that have time variability dif- 
ferent from that of 60 Co. In particular, the prompt 
background lines of 69 Ge at 1337 and 1347 keV show 
up as strong negative residual lines in the spec- 
trum. Note that also at 1173 keV residual counts re- 
main, while at 1332 keV the spectrum appears well- 
represented (insignificant residuals). 



4. SPECTRAL ANALYSIS 



Straight background spectra subtraction does not 
make use of the coded- mask and field-of-view proper- 
ties of SPI; on the other hand, image deconvolution is 
feasible only for strong signals. As we expect a weak 
or no signal, we derive spectra by fitting the ampli- 
tude of a given skymap of (expected) intensity dis- 
tribution, per 1 keV wide energy bin. As skymap we 
used the Dirbe 240 fim intensity distribution, since 
this map has been shown to provide a reasonable 
tracer of galactic 1809 keV emission, and hence is 
probably also a r easonable first-ord e r trac er of the 
60 Fe distribution (|Knodlseder et all 11999ft . Fig. H 
shows the resulting spectra. 

As in the residual spectrum (cf. Fig. 0J), a signifi- 
cant line feature appears at 1173 keV. The formal 
flux in this feature amounts to (9.3 ± 2.6) x 10~ 5 
ph cm _2 s _1 , which corresponds to ~ 30% of the 
1809 keV line flux. This is thr ee times as mu ch as 
the flux reported by RHESSI ijSmitbt l2004a|) 1 and 
also exceeds prev i ously reported upper flux limits 
ipiehl fc Timmesl 11998ft . At 1332 keV no line fea- 
ture appears. The formal 1332 keV line flux amounts 
to (-0.2 ± 2.4) x 10~ 5 ph cm-V 1 , 



5. DISCUSSION AND CONCLUSIONS 



If the feature at 1173 keV were a real celestial signal, 
we would have expected to find the same flux values 
for both lines. In addition, the formal flux in this 
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Figure 3. Residual SPI energy spectrum after subtraction of the background model (cf. Fig. QJ). 
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Figure 4- SPI spectrum of the galactic centre region (GCDE1 and 2), assuming an 60 Fe line intensity distri- 
bution that follows the Dirbe 240 /im infrared emission as a tracer of the massive star population. The ordinate 
gives the 60 Fe line flux in units of 10 -4 ph crre - 2 s — 1 keV -1 for the central steradian of our Galaxy. The residual 
of the instrumental background line at 1337 keV from 69 Ge has been modelled by fitting a Gaussian function to 
the spectrum and by subtracting this function from the data points. 



features exceeds the expectations based on measure- 



ments by other instruments. Therefore we conclude 



that the spurious line residual at 1173 keV, inconsis- 
tent with 60 Fe due to the absence of a corresponding 
signal at 1332 keV, should be used to determine our 
systematic uncertainty at this point. From this, we 
derive an upper flux limit in each of the 60 Fe lines 
of about 1CP 4 ph cm~ 2 s _1 . In comparison, the sta- 
tistical uncertainty of ~ 2.5 x 10~ 5 ph cm _2 s~ 1 is 
small. 

Our current data analysis efforts are concentrated on 
reducing the systematic uncertainties and improv- 
ing the subtraction of the instrumental background 
lines. Once the systematic uncertainties have been 
reduced to a reasonable level (i.e. below the sta- 
tistical uncertainties of the data) we can expect 2a 
upper flux limits of the order of 5 x 10~ 5 ph cm _2 s _1 
in each of the 60 Fe lines for the first year of GCDE 
data. Combining both lines would reduce the up- 
per flux limit by roughly \/2, leading to a level of 
3.5 x 10~ 5 ph cm -2 s -1 . Exploiting also the multiple- 
event data (which obey about the same signal-to- 
noise ratio as the single-event data) will give an- 
other factor of y/2, providing an upper flux limit of 
2.5 x 1CT 5 ph cm" V 1 . This is about 6 - 8% of the 
1809 keV line flux for each of the 2 lines (dependent 
on which flux level is assumed for the 1809 keV line 
intensity) and would present the most stringent up- 
per flux limits available today. Adding more data to 
our still rather small dataset (only 3.8 Ms of exposure 
time have been exploited so far, yet much more data 
have already been collected along the galactic plane 
by INTEGRAL) will further reduce the sensitivity 
limit. 

On the other hand, as stated in section ^ cosmic 
abundance arguments suggest a 60 Fe line flux of 
~ 3 x 10~ 5 ph cm _2 s _1 in each of the lines. This 
is comparable to the 2a upper limit that is achiev- 
able by SPI using the first year's GCDE data and is 
close to the value suggested by RHESSI. Adding all 
INTEGRAL exposure of the central galactic regions 
from the first mission year will probably bring the 
expected 60 Fe line intensity near to a firm detection. 
Thus, it is highly likely that the spectrometer SPI on 
INTEGRAL will provide in the near future serious 
constraints on 60 Fe nucleosynthesis in our Galaxy. 
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